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PREFACE. 

TI1HE lectures contained in tliis booklet were all delivered 
practically extempore. That they s<^ 1 the light of day, is 
chiefly due to the kindness of Mr. Corakh Prasad, M.Sc , and 
Mr. iBadrinath Prasad, M.Sc., Assistant Professors of Mathematics 
in the Benares Hindu University, of whom the former took down 
notes of the last three lectures and the latter of the first three. 
The copy for the press was prepared by me from these notes 
with the help of Mr. Mewa Ram Saksena, M.A., H.Sc., of the 
Shyamsuiidar Memorial Intermediate College. Chandausi. Dr. 
Bibhutibhusan Datta, D.Se., Acting Head of the Department of 
Applied Mathematics in the University College of Science, Cal- 
cutta, kindly helped me in correcting the proof sheets. I take 
this opportunity to express my thanks to all these friends and 
former pupils. 

OANESH PRASAD. 

April. 1924. 
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FIRST LECTURE. 

[Delivered on the I Oth January, 1921.] 

INTRODUCTION. 

The chief reason for my choosing " The place of partial 
differential equations in Mathematical Physics" as the subject for 
these lectures is my wish to inspire in my audience a love for 
Mathematics. Before entering into details, however, I shall give a 
brief historical account of the application of Mathematics to natural 
phenomena. 

Mathematicians devote themselves to their subject chiefly with 
one of the two objects : either to study natural phenomena in con- 
nection with matter, space or time ; 01 for purely aesthetic pleasure. 
A good illustration of the latter is afforded by their attempt to find 
all the prime numbers less than a certain given number. In ancient 
times people studied Mathematics mainly because of some mystery 
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which they attached to it and not for its application in studying 
natural phenomena ; in this connection, we may note, for example, 
the doctrines of Pythagoras or the Sulvasutras of Apastamba. 
Among the earliest writers who devoted themselves to Mathematics, 
were Archytas of Tarentum (428-327 B.C.), Aristotle (3 84-3 2 2 B.C.) 
arid Archimedes (287-212 B.C.). Archytas applied his knowledge 
of Mathematics in building up the science of Mechanics. After 
Archytas, Aristotle did something in this respect. But he had 
strange ideas about some elementary truths ; for example, he 
supposed that for a falling body the rapidity of its motion depend- 
ed upon its weight, according to him a heavy body falling more 
rapidly than a lighter one. Another erroneous idea which Aristotle 
had was that if a stone were thrown up at a certain angle, it 
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would go first upwards in a straight line at that angle and then 
fall vertically downwards, the path of the projectile consisting of 
two straight lines, as in the accompanying figure. Archimedes 
discovered the principle of the lever and the conditions of the 
floatation of a solid in a liquid. Galilei (1564-1642) was the first 
to acquire clear ideas about the first two laws of motion. He also 
stated and proved that the path of a projectile was not a pair of 
straight lines but a parabola. After Galilei, Mechanics became 
gradually a science. A large number of results were discovered 
between 1642 and 1700. 

In the beginning of the 18th century a new epoch was created ; 
for, people began to think of natural phenomena other than 
ordinary motions, such as for instance the state of a hot body. 
This gave rise to the use of partial differential equations. Let there 
be a solid body which conducts heat. The amount of heat passing 




FIG. 2. 



across a sheet at P perpendicular to the ic-axis is required to be 
found out first. The amount of heat conducted will be propor- 
tional to the gradient of the temperature and is equal to 



which is the first germ of truth in the Analytical Theory of Heat 
and was discovered about the year 1747; the partial differential 
equation governing the state of a hot body was discovered much 
later. * Another important discovery was made in connection with 
the study of sound. People had already acquired some idea of 
harmony, and the formula of Taylor, 



viz. 
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was roughly known to Aristotle and the ancient Riahis of India. 
In 1636 Mersenne gave the above formula as the result of his 
experiments. But all such results were at first more or less of an 
Bmpirical nature. The equation arising in the problem of the 
vibrating string, viz. 

ty-T^y 

di* ~ P ft? 

was given by D'Alembert as lato as in 1747 and was the first partial 
differential equation to be carefully studied. It was solved in 

T 

two ways. Putting = a 2 , 

y = f(x + at) + F(x - at) , 

was the first solution given by Euler and D' Alembert. Another, 
which was given by Bernoulli, is 

* . . UirX YlTTdt 

y = 5 A n sin cos ~ . 
! / / 

]*n appearance, Bernoulli's solution is periodic whereas the 
former solution is non-periodic. Therefore the two solutions 
formed the subject of a most interesting controversy. Fourier 
showed that Bernoulli's solution was also right and that, contrary 
to the views of Euler, it can be used for the case of an arbitrary 
displacement. 

It is the deliberate opinion of modern authorities that before 
the use of partial differential equations, there was no sound basis 
upon which the science of Physics could rest. Soon after the time 
when partial differential equations came to be applied to the study 
of natural phenomena, Fourier made the following statement 
which will remain true for ever : 

" The analytical equations, unknown to the ancient geometers, 
which Descartes was the first to introduce into the study of curves 
and surfaces, are not restricted to the properties of figures, and to 
those properties which are the object of rational Mechanics ; they 
extend to all general phenomena. There cannot be a language 
more universal and more simple, more free from errors and from 
obscurities, that is to say more worthy to express the invariable 
relations of natural things. 

" Considered from this point of view, mathematical analysis is 
as extensive as nature itself; it defines all perceptible relations, 
measures times, spaces, forces, temperatures ; this difficult science 
is formed slowly, but it preserves every principle which it has once 
acquired ; it grows and strengthens itself incessantly in the midst 
of the many variations and errors of the human mind. 

" Its chief attribute is clearness ; it has no marks to express 
confused notions. It brines together phenomena the most diverse, 
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and discovers the hidden analogies which unite them. If matter 
escapes us, as that of air and light, by its extreme tenuity, if bodies 
are placed far from us in the immensity of space, if man wishes to 
know the aspect of the heavens at successive epochs separated by 
a great number of centuries, if the actions of gravity and heat are 
exerted in the interior of thr earth at depths which will be always 
inaccessible, mathematical analysis can yet lay hold of the laws of 
the phenomena. It makes them present and measurable, and 
seems to be a faculty of the human mind destined to supplement 
the shortness of life and the imperfection of the senses ; and what 
is still more remarkable, it follows the same course in the study of 
all phenomena ; it interprets them by the same language, as if to 
attest the unity and simplicity of the plan of the Universe, and to 
mako still more evident that unchangeable order which presides 
over all natural causes." 

We have the following important partial differential equations 
that occur in the study of natural phenomena : 



.. 

( h) + - + - = tira. 
' <9:r* fly 2 flz' ' 

' f)V fl'F 



As an illustration of the development of the rigorous method 
for the study of natural phenomena, I will sketch an outline of 
the history of the proof of Poisson's equation, 



The proof of Poisson's equation, usually given in ordinary 
text-books, is not quite correct. The proof is based on enclosing 
the point (x, ?/, z) in a small sphere where the density t> is taken to 
be constant. F= F, + V where T 7 2 is the potential due to the small 
sphere and F, is the potential due to the rest of the attracting 
body. v*F,=0 and thus v*F = r 2 F i . Now if you postulate that p 
is constant everywhere in the small sphere, then of course 



FIRST LECTURE. 5 

proof can hold true. But in case p IB variable, this proof will no 
longer hold. Green took up this question in 1828 and said that p 
is not constant but its variation is very small in the sphere and is 
therefore negligible. The first rigorous proof was published in 
140 by Gauss who took p and its first partial differential coeffi- 
cients to be continuous. Now suppose F, the potential at P(x, y, z), 
is due to a solid for which / at any point Q is 

cos*0 



where PQ = r, the radius vector and is the angle made by the 
radius vector with the axis of x. Tn this case Poisson's equation 
is not valid at P ; for, 

<9 2 F d*y 



has no meaning. I invite your attention to Professor Petrini's 
generalization of Poisson's equation, viz. 

+ k,z ay 



dy 



Sz dz 



where A, k, and I tend to zero without the ratio of any two of 
them tending to zero or infinity. There are various advantages of 
this generalization, but in my paper l " On the failure of Poisson's 
equation and of Petrini's generalization/' I indicated a case in 
which this generalization fails. Take for example 



then there is no limit in this case except for very special values of 
A, k, 1 9 and the generalization as well as Poisson's equation fails. 

The application of the modern Theory of Functions of a real 
variable makes the use of partial differential equations more logical. 
Also the application ot the Theory of Integral Equations is very 
useful. A typical integral equation is 



Bulletin of the Calcutta Mathematical Society, Vol. 8, 1918. 
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where <f> is the unknown function to be determined. The first im- 
portant paper which laid the foundation of the Theory of Integral 
Equations as a systematic branch of Mathematics was published 
by Professor Ivar Fredholm in 1900 and was later published in the 
Acta Maihematica in an extended form in 1903. Integro-differential 
equations of a somewhat complicated type are useful in the " Me- 
chanics of Heredity," i.e. those branches of Mechanics in which the 
.state at a particular point of time is dependent not only on the 
immediate past but also on the whole period antecedent to that 
point of time; as for example in the phenomena of Hysteresis. 
Intogro-differential equations are equations which involve not only 
the unknown functions under signs of integration but also the 
unknown functions themselves and their derivates. A paper on 
the subject of integro-differential equations by Professor Volterra 
appeared in the Acta Mathematica, Vol. 35. The subject of integral 
equations is of essential importance in studying the phenomena of 




Fio, 



radiation. The law of Kirchoff, - = constant, cannot be proved 

a 

without the help of integral equations. The first rigorous proof 
was given by Professor David Hilbert in 1913. The kinetic theory 
of gases cannot be rigorously studied without the help of integral 
equations (see Hilbert's paper in the Mathematische Annalen, Vol. 
72, 1912). 

There can be no better instrument for studying natural 
phenomena than integral equations, partial and ordinary, and 
integro-differential equations. The truth of this remark can be 
easily seen by considering the problem of Dirichlet : To find a 
solution of v*F = inside a given region (Fig. 3) when, on the 
surface bounding the region, F=/ (x, y, z), a prescribed function. 
By the help of the integral equation 
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the quantity p can be found as the strength of a double layer whose 
potential is the required solution of the problem. Integral equa- 
tions are also most serviceable in studying the problem of the 
conduction of heat in a body. The well-known partial differential 
equation for the linear conduction of heat is 

dV 



But this equation is invalid for some very simple cases and has to 
make way for a partial integral equation. (See G, Prasad's memoir 
" Constitution of matter and analytical theories of heat " in the 
Abhandlungen der Koniglichen Gesellschaft der Wissenschaften zu 
Gottingen, 1903). 



SECOND LECTURE. 

[Delivered on the llth January, 1921.] 

THE PARTIAL DIFFERENTIAL EQUATIONS OF VIBRATORY 

PHENOMENA. 

The most important partial differential equation met with in 
the treatment of ordinary vibratory phenomena is of the form 

r a*F a*r a 2 



Vibratory phenomena of electrodynamio origin are governed 
by the equations of electrodynamics : 



By dp 

4:itU = ~' --- r- , 

a/ a* 

da dy 

47ri> = - -- 3r , 

a? a# 
a/? a a 

4:7TW = -T~ 7T- 



dy 



dc Y dX 
and ~^=T"~5~" 

a^ dj- a^/ 

where the components of magnetic force at the point (x, y, z) are 
denoted by <*, ft, y, the components of magnetic induction by a, b, r, 
the components of electric force by X, Y, Z, and the components 
of total current strength by u, v, w. These six equations represent, 
in symbols, the laws of Faraday, and their solution reduces to the 
solution of a single partial differential equation of the type 

A 4 

* 

In the study of many kinds of elastic vibrations, we come 
across a partial differential equation of a similar type. In the case 
of the vibrations of an elastic plate, the motion is given by solving 
an equation of the form 
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"<- 

<9 >z a* & 

where, as usual, y* stands for + + -r-j . 

I proceed to give you two examples to show how partial 
differential equations are utilized to determine the vibrations in 
concrete cases : 

First consider the motion of air in a rigid spherical envelope. 
Transforming the equation 



into polar co-ordinates, we have 



d ' l 

8m 



- - 

df ~~ r* 1 flrV flr / sin 30V 90 / sin*0 cty* 

Putting F=r/ (r, 0, <p) . cos kct, we find that the equation is 
satisfied by 

J *<*(&) 



where S n (0, </>) is any surface harmonic of degree n. Certain bound- 
ary conditions are to be satisfied in every case. Here 



a being the radius of the envelope. This gives the suitable values 
of k and thus we have an infinite number of solutions each satis- 
fying the differential equation and the condition at the boundary. 
The sum of such solutions 

WV) 

F= 2 - ^~S A (0, <J>) cos k ct 



A 



must also satisfy the differential equation and the boundary con- 
dition and gives the motion, the constants in $ W|A being HO chosen 
as to give the prescribed initial state when t is equated to zero. 

As a second example, consider the vibrations of a circular 
membrane. 



is the equation of motion. On transforming it into polar co-or- 
dinates, we have 



10 DIFFERENTIAL EQUATIONS IN MATHEMATICAL PHYSIOS, 

Proceeding as in the first example, 

I dU 



The typical solution of this equation is 

7 = J n (kr) (A n cos n + B n sin n 9). 

Since at the boundary r = a, there is no displacement, U is 
zero there and therefore k is a root of the equation 



The motion is given by 

F= 2 Jn(kfj) cos & A ct (A n A cos nO + B n ^ sin w#) , 

n, \ 

the constants being determined by the prescribed initial state. 

Now I want to give you some criticisms, from the point of 
view of the theory of functions of real variables, on the usual 
method adopted for the study of vibratory phenomena. For the 
sake of simplicity and fixity of ideas, I will consider the simplest 
vibratory phenomenon, viz. the vibrations of a string. 

The partial differential equation of the motion is 



df djf ' 

T 

where c* stands for and may be taken to be unity if the units 
p 

are suitably chosen. Let the string be fixed at its extremities^ 
#=0, #=7r and let it be set in motion by a displacement ,y=/ (x). 
Then the motion at any subsequent time t is given by 



y = 2 ^4 w sin nx cos nt , 

=] 

the constants being chosen in such a way as to satisfy the initial 

condition . Therefore 

f 
f(x) = 2 .4 M sin nx 



and ^, = -1 f(x) sin nxdx. 



The flaws in the above method for determining the motion of 
the vibrating string were pointed out by me in a paper l " On the 



See Vol. 6 of the 4 < Bulletin of the Calcutta Mathematical Society," 1916. 
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vibrating string with an infinite number of edges." The following 
are the two flaws : (1) We always assume the existence of 

dV * 
*? and A? 

rft every point. (2) It is taken for granted that if each term of an 
infinite series satisfies the differential equation then the sum of the 
series will also satisfy the equation. 

Both the assumptions are unjustifiable. For, a string may 
have an edge so that at that point 

, what to speak of , 
ox ox* 

will be non-existent ; also what is true of a finite number of terms 
may not be true of an infinite number of terms. These objections 
I proceed to elucidate at some length. 

Consider the string which is under constant tension T. Ac- 
cording to Newton's second law, study the change of momentum 
of a small element of the string in the period of time 



Let the #-co-ordiriates of the ends of the element be x { h> o^ 
Then the change in its momentum is 



It must be balanced by the change in the force which is 

'l+T 

' //)*/ v ,/J,.iv ^ 

!. (1) 



Only the assertion of this equality gives the true equation of the 
phenomenon as an integro-differential equation. 

What is usually done is to consider 2 h as getting smaller and 
smaller and also 2 T getting smaller and smaller. Then the appli- 
cation of Rolle's theorem shows that the integro-differential equa- 
tion is equivalent to 



the 0's are all proper fractions. 
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The conclusion, that 



is not always right ; if there is discontinuity then the above equa- 
tion may be invalid. 

Now I proceed to elucidate the second flaw by taking the 
initial state of the string to be f/ = /(.r), where 



- , and / (x)=w-x for ^ 



The ordinary procedure gives the displacement y at any time as 

4 5j _] sin (2w. 1)# cos (2n-l)J 

" ' ' 



the assumption being that 



sin (2w-l);c cos (2ft- I)/, 



and ^-?=~~ 2 (-I)"' 1 sin (2?i~ l).r cos (2n-l) 

oa?' TT w= i 

and therefore the equation 



is satisfied. But the assumptions are unjustified because the series , 
to which 

7TY as well as - 
is equated has no meaning, it being not convergent. Thus both 

^ and ^l 

have to be investigated by a different procedure and their equality 
established by a different method. 

How to study then the phenomenon of the vibrating string ? 
This question was answered correctly by the late Dr. A. Harnack * 
in 1887. If 



l See his paper, "Ueber die mit Ecken behafteten Scbwingungen gesoannter 
Sait6n" (Mathematische Annalen, Bd. 29). 
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are finite and continuous then it was proved by Harnack that the 
most general motion of the string is given by 

oo 

y=z % A n sin nx cos nt. 

n = l 

Proof : 

Since ~ is existent and finite, according to Fourier's theorem 
dx 

y (or, t) must be equal to 

00 

S a w (t) sin nx, 

w = l 

7T 

2 
where M =- I y (x, t) sin w# dx. 



[Here it may be noted that every continuous function is not 
expansible in Fourier's series. Thip was first shown by Du Bois 
Reym<jnd who gave an example of a continuous function the 
Fourier's series corresponding to which is not convergent.] 

Integrating by parts 



7T 7T 

r rj" 2 r 

I ^ sin nx dx = w* \ y si 

.1 dx y 



sin nx dx^ n % - a n (t). 



f 

Since ^ is continuous, we have 

dr 



a f . 

- I v s 
81) 



. , . 

sm ?i.r dr= \ ~ BUI nx dx, 



(I 

which equality is true because ~ is continuous. Therefore similarly 

nt 



a * (' 

"fa \ y s ' 


( 

and this latter integral equals 



j 

nx dx= \ ;-j- sin 



-r-r sin 
fix* 



, 

because 



Hence *- a w - n^a,, = 0. 
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Therefore o>n A n cos nt + B n sin nt. 

If the string is simply displaced initially, then from 

00 

y 2 (A n cos ?i + B n sin n/) sin no: 

tt = 

we find that the B's must be all zero. 
Therefore the general solution is 

X 

y= 2 ^ w cos n< sin n,r. 
i 

1 proceed to consider the case when there is an edge at a point 
P on the string so that at that point ~- is non-existent and conse- 

quently r-~ is also non-existent there. 

^ J ax* 

Denote the ordinates by ?/,, y. 2 in the portion APB and by 
y\, y\ in the portion AP'B. The usual phenomenon of wave mo- 
tion is that every wave propagates itself. Consider the small ele- 
ment of the string round P during a small period of time r. This 




N 



B 



FIG. 4. 



. dc 



element is * where c is the abscissa of P ; the element undergoes 
dt 

a change 

dc c dz/j dt/ 2 i 

P Jt T ( ~di~Hit \ 
in momentum and the change is brought about by the impulse 



acting for the period of time 
being suitably chosen, at P 



T 
Therefore, since - = 1 , the units 



_ 

dt dt to 



_ 

dx dx 
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As the string is continuous y { (c, t) = j/ 2 (c, t) ; therefore 

y { dc dy\ 9t/ 2 dc dy% 
Ic * "^*"dF = !c~* Te^TST' 

This gives &.&*(^.^ (3) 

In order that the two equations (2) and (3) be consistent, 

at 

must be equal to 1. Therefore the equation, which must be sa- 
tisfied at the edge in the place of the partial differential equation, 
which has no meaning there, is 



d l\ . ^ 

dt/c-o Vdt/"" VfljJ/e + o \dx/e-o ' 



l d l 

\dt 



The question is now whether in the present case the general 
solution holds ; it will bo easily seen that the required condition 
at the edge is satisfied by the solution and therefore the answer 
to the question is in the affirmative. 

Cases of strings with three or four or any finite number of 
edges have been considered by Harnack. But the case of an in- 
finite number of edges was treated for the first time by me in the 
paper already referred to. As a simple example, it is proved there 
that 

2 On 

y~- 2 -~r cos nt sin nx 

TT n^tf 

gives the motion of the string, vibrating with the infinite number 
nf edges, the abscissae of which initially formed the set 



G n being equal to 

n M-l 

, ^i" 1 or, 1 . nw ' ! nv 

(- 1 ) 2 7 8in 7T or(-l) 2 7<*>s T , 

i=l & & 4=1 *J ^ 

according as n is even or odd. 

I will conclude this lecture by showing how the theory of 
integral equations may be applied in the study of the problem of 
the motion of a vibrating string. 

Let us consider the case of a string of variable density. Then 
in the equation 



P IB not constant. Lord Rayleigh could not integrate this equation 
except for very special types of P. 
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The general solution can be obtained by the help of an Integra, 
equation. For, put 

2/ = cos nt . U(x). 

Then the equation (4) becomes 

-^+y,rt*) */<>, 
(jc j. 

and it can be shown without difficulty that this last equation to- 
gether with the boundary conditions 



is equivalent to the integral equation 



where G (f. x) - f_ f or 
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f Delivered on the 12th. January, 1921.] 

THE PARTIAL DIFFERENTIAL EQUATIONS OF THE ANA- 
LYTICAL THEORY OF THE CONDUCTION OF HEAT. 

After the discovery of the experimental foundations of the 
elementary laws of conduction and radiation by Newton, Lambert 
and Biot, a good deal of time elapsed before the discovery by 
Fourier of the following differentia] equation as representing the 
phenomenon of the conduction of heat in an isotropic solid : - 

<?Z_ 

17 ~~ c 

where *c 2 is the coefficient of diffusivity. A number of special 
problems were solved by Fourier and his contemporary Poisson. 
For example, take the problem of the cooling of a solid sphere, 
whose surface is maintained at a constant temperature, say zero : 
the initial state being given by 

V(x, y, z,0)=f(x,y, z) . 

The particular case in which f(x, y, z) is a function of r alone 
was solved by Fourier ; the general case was considered by 
Poisson. To determine the state inside the sphere at time t, we 
begin by transforming the partial differential equation into polar 
co-ordinates. Thus 



dv - *[~ l i d f * dv 
~dt- C LrM STrV l^ 
Putting 



we have 
d 



of which the typical solution is 

* fe") o (6 

o n (0, 



where S n (0, </>) is a spherical surface harmonic of degree n. At the 
boundary r=a, F is zero ; therefore k must satisfy 

/. + * (fo)0. 
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Let the values of k be A/ p & 2 , k.^ etc. Then the required solu- 
tion of the problem is given by 



the summation extending over all values of n from upwards and 
over all values of m from 1 upwards and the constants in S n , m being 
so chosen as to satisfy the condition relating to the initial state. 

As another example, consider the state of heat of an ellipsoid 
of three unequal axes. Larne in 1837, Heine in 1839 and Mathieu 
in 1870 made important contributions to the advancement of the 
knowledge of the subject. Transforming the partial differential 
equation V 2 F = into ellipsoidal co-ordinates, the equation 
reduces to 



where Re stands for \S(a* + 6) (6* + 6) (c*~+J) . Taking F td be of 
the form E(h) E(/u) E(v) where E(0) is a function to be determined, 
the equation will be satisfied if 



where A and B are constants. The function B(0) is called Lame's 
function and the product E(\) E(f<) E(v) is called an ellipsoidal har- 
monic. The stationary state of heat of* the ellipsoid can be 
expressed in terms of a series of such products. As regards thd 
non-stationary state of heat, the partial differential equation of 
the phenomenon is 

9 _Z=c*v*F 

dt 

and this reduces to the form 



*<>-(*)'}- 

where 



Mathieu made an unsuccessful attempt to solve the above 
equation. For the special case of an ellipsoid of revolution, Professor 
C. Niven's results are of some interest. The most successful 
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attempt to solve the problem of the non-stationary state that has 
been made up to now is that of Mr. Bibhutibhuson Datta, who 
obtained, for the general case of the ellipsoid of three unequal 
axes as well as for the particular case of the ellipsoid of revolution, 
approximate results of very great interest and importance, one of 
tl\p results being the correction of a mistake of Prof. Niven. (See 
the Bulletin of the Calcutta Mathematical Society, Vol. 8, and the 
American Journal of Mathematics, Vol. 41). 

T propose now to submit the use of partial differential equa- 
tions to a criticism from the point of view of one who wishes to 
be strictly rigorous. For the sake of simplicity and fixity of ideas, 
I will consider the conduction of heat in an infinite slab bounded 
by the faces x= - TT, # = 71-, which are impermeable to heat, and T 
will take the initial temperature of the slab to be f(x) where f(x) 
is an even function of x. 

The partial differential equation of the phenomenon is 



- 

dt ~~ G dx 1 ' 

where c* stands for the coefficient of diflfusivity and may be 
taken to be unity if the units are suitably chosen. Also because 
of the impermeability of the bounding faces, #= v, x = ?r, we 

dV 

have = for x = v as well as for # = - T. Thus the state of the 
dx 

slab is given by 



'where 



00 rfit 

V = JU; M + 2 a n e ' cos nx, />^0, T 



9 
7- w = ~ \ f( x ) cos nx ( l' x * 



The above-mentioned treatment of the problem is the ortho- 
dox one and will be found in every text-book. In Fourier's 
original papers the statement of the conditions is the same as 
above, with this difference that it is clearly mentioned that 



for every positive value of t including zero and for every value of 
x in the interval ( - w, w). 

In order to establish the defective character of the orthodox 
treatment, I take the simple case in which f(x) is given by the 
accompanying graph. 
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At there is no differential coefficient and still by the ortho- 
dox method 




-tT 




Fio. 5, 



Km 



~" 



If we differentiate the abqve 
Aeries term by term the series 
so obtained is uniformly con- 
vergent as long as t > and 
hence we can differentiate in 
that manner for such values 
of t. Thus for t > 0, 



dt , 

and, making t tend to 0, 

= oo . 



Now according to the theory of functions of real variables, if 
lim dV(0,1) 



is existent, it muet be equal to 

lim F((U)-F(0,0) 

Therefore, at the time = () the value of -=- at is infinite 

ot i 

and also at that point, there being no differential coefficient of F 
with respect to x, v-^-is non-existent ; consequently for # = 0, =0, 
the equation 



has no meaning and the orthodox treatment fails. 

How to study then the phenomenon of the conduction of 
heat? 

This question was answered correctly for the first time in my 
memoir, " Constitution of matter and analytical theories of heat ", 
published by the Royal Society of Sciences of Gottingen in 1903. 
(See Abhandlungen der Koniglichen Gesellschaft der Wissenechaften 
zu Gottingen, Bd. 2, No. I). In the theory, as remodelled by me, 
the equation of the phenomenon is shown to be not a partial 
differential equation but a partial integral equation, and attention 
is drawn to the simple principle that it is physical continuity 
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which is necessary and not the continuity of the form of an ana- 
lytical expression. I proceed to elucidate these features of the 
theory. 

First, as regards the equation of the phenomenon, this 
follows from the application of the principle of the conservation 
'of energy in accordance with which the amount of heat passing 
into a volume must be equal to the amount of heat absorbed to 
raise the temperature of the mass of the volume. In the form 
of symbols, this condition may be written as : 



(A) 



\ d ^l^ldt'=\ ^ 



T(x 9 t) denoting the temperature. It can be easily shown that 
this integral equation is equivalent to the differential equation 



^ 
ftt ~fe* 

dT d*T 

if _ and - are finite and continuous. For, by the mean value 
dt ox* 

theorem, 



i 

) 



and ^ 2 being proper fractions. Similarly 



From (A), the two integrals given above are equal; consequently, 
proceeding to the limit by making 2 - t l and x. l - x { tend to zero, 
we have 



for =aj|, tsst l9 if these differential coefficients are finite and 
continuous. 

Secondly, T(x, 0) need not be derivable from the analytical 
expression, which equals T(x, t) when t> 0, by putting ^=0 in that 
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expression. Thus, in the remodelled theory the convergence of 
the Fourier's series 

00 

|<7 + 2 a n COB nx 

i 

has no place. On the other hand the continuity of T(x, t), i.fe., 
the equality of 

lim V(x. t) to /(a?), 
= 4-0 

follows at once from the continuity of f(x) ; and the only restric- 
tion imposed on f(x), in order that V(x, t) should be the solution, 
is the following : 

(B) For each value of #, within the interval (-TT, *), the upper 
<9F 



limit of 



in the neighbourhood of = 0is finite. 



Sx 

SV 
The necessity for the finiteness of T~ is obvious ; for, since 

J ox 

dV 

- k is the amount of heat which flows across an unit area of the 
de- 

plane at x in unit time, if this amount becomes infinite it will 
be absurd. Now for >(), the series for V(x, t) being uniformly 
convergent and the series obtained by differentiating it being also 

dV 
uniformly convergent, ^- is finite. Hence we have to consider only 

the neighbourhood of = 0. 

The equation (A) is satisfied by 

T(x> t)=V(x, t),t>0 
as long as f(x) is finite and integrable. For, for />0, the series 

for "2 and V can be integrated term by term, the series concern- 



As regards (B), it can be satisfied even if -- makes an infinite 



ed being all uniformly convergent. 
As regards (B), it can be satis! 
number of oscillations as t tends to ; for example, if 

/(*) 



)=! cos (log ^j dx , 




dV 
gr makes an infinite number of oscillations with indefinitely 

increasing amplitudes as t tends to zero. Thus theoretically the 
slab may, at certain points, become hotter and cooler an infinite 
number of times. 
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I conclude this lecture by pointing out that a number of in- 
teresting problems in the theory of the conduction of heat can be 
best solved by the use of the theory of linear integral equations. 
For example, one of the earliest applications of that theory is by 
Professor Picard, to the problem of the cooling of a ring. The 
partial differential equation of the phenomenon is 



This differentia] equation can be solved in terms of the solu- 
tions of an ordinary differential equation of the second order which, 
in turn, is equivalent to an integral equation of the type 



U(i) - F(() + XI U(x) A Lx) G(x, Qdx , 
where TJ is the unknown function. 



FOURTH LECTURE. 

[Delivered on the 3rd February, 1921.] 

THE PARTIAL DIFFERENTIAL EQUATIONS OF THE THEORY 

OF GRAVITATIONAL ATTRACTIONS AND THOSE OF 

ELECTROSTATICS AND MAGNETOSTATICS. 

In the three lectures which I delivered last month, I first gave 
a historical review of the present theory of the partial differential 
equations of Mathematical Physics together with some criticisms 
and then dealt in succession with the equations of vibratory 
phenomena and of the analytical theory of the conduction of heat. 

I wish to preface the three lectures which T intend to deliver 
in this month by explaining to you the two kinds of Mathematical 
Physics. This cannot be better done than in the following ^ords 
addressed by Professor Vito Vol terra, of Rome at the inauguration 
of the Rice Institute of America in 1915 : " There are two kinds 
of Mathematical Physics. Through ancient habit we regard them 
as belonging to a single branch and generally teach them in the 
same courses, but their natures are quite different. In most cases 
the people who are greatly interested in one despise somewhat the 
other. The first kind consists in a difficult and subtle analysis 
connected with physical questions. Its scope is to solve iri a 
complete and exact manner the problems which it presents to us. 
It endeavours also to demonstrate by rigorous methods statements 
which are fundamental, from mathematical and logical points of 
view. 

" I believe 1 do not err when T say that many physicists look 
upon this mathematical flora as a collection of parasitic plants 
grown to the great tree of natural philosophy. But this disdain 
is not justified. In the evolution of Mathematical Physics these 
researches are to play an ever increasing part 

" The other kind of Mathematical Physics has a less analytical 
character, but forms a subject insuperable from any consideration 
of phenomena. We could expect no progress in their study 
without the aid which this brings to them. Could any one imagine 
the electro-magnetic theory of light, the experiments of Hertz and 
wireless telegraphy, without the mathematical analysis of Maxwell 
which was responsible for their birth ? " 

To-day I shall confine myself to the first kind of Mathematical 
Physics and the lecture will be mainly destructive in character. 
But to-morrow I shall speak on the constructive side of things, 
when I shall take up the second kind of Mathematical Physics. 
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The following are the most important equations of gravita- 
tional attractions, Electrostatics and Magnetostatics with the 
respective years of their discovery : 



(1) 0T + a"T + Tp- * (Laplace's equation, 178 1) 

(2) TT + - = - 47m, (Equation of Coulomb and 
dn l dn. z ^ * - 

Poisson, 1811) 
being the surface density. 



(3) P + - -1- = -4*7*, (Poisson's equation, 1813) 

p being the volume density. 

(4) W^- tt 7 _ = -47m, (Helmholtz's equation, 1853) 

a being the strength. 

[T^he exact form of the equation (2) was given by Poisson in 
1811. The last equation gives the potential difference at the 
surface of a double layer of strength -.] 

Let us take up Poisson's equation 



3-y ~7 FT" ~ ~ ' 
dor at/ ^ a^ z 

The proof as given in text-books postulates /> to be a constant 
throughout a sufficiently small sphere round (#, y, z). Let u and 
f7 denote respectively the potential due to the small sphere 
round (x, y, z) and that due to the remaining mass. Then 



By Laplace's equation V * U = 
stnd consequently r a F= r 

Taking /o to be constant 



at a point (f, >;, ) in the sphere distant r from (^; ?/, z) 3 a being the 
radius of the sphere. Differentiating two times with respect to 
y 9 f we have at the centre 

?*=-?!! (2 + 2 + 2) =-4^. 

o 

In 1828, Green gave a proof more or less like the above. The 
first really rigorous proof was given by Gauss in 1840. Gauss's 
proof, which I will reproduce presently, is based on the supposition 
that p is finite and continuous together with its three partial 
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differential coefficients of the first order. In 1882, Professor Holder 
succeeded in establishing Poisson's equation on the supposition 
that 



(, 77, )-p(x, y, z) 



where i* and A are positive constants and r is the distance between 
(, v> ) and (a;, y, z). The next investigator to impose still less 
restrictions on p was G. Morera who proved in 1887 that, if 






is finite, then Poisson's equation is valid. In 1891, Kronecker 
attempted to establish the equation with losser restrictions than 
the preceding investigators. From 1891 to 1900 practically nothing 
was done to improve the proof of Poisson's equation. In 1900 and 
the next ten years, Dr. Henrik Petrini tried to make the proof as 
little restrictive as possible. Some of my papers deal with the 
criticism of Petrini's work. Petrini 's chief result is that if 

2vr 1 a 



im I d<t> \ (3//-l)d/. \ 
= <>J J J r 



lim 
h 

-J 



u 

exists and is finite, then ^ exists and is finite, p standing for 

c/x 

cos where is the angle between the radius vector and the axis 

rva r\r> 

of x - with similar conditions for -r> and TT-T ; the additional con- 

dy l dz 2 

dition for the validity of Poisson's equation being the continuity 
of /> which condition, however, has not been sufficiently emphasized 
by Petrini. But cases may occur in which p is discontinuous. 
For example, I pointed out in a paper published in the Philo- 
sophical Magazine in 1917 that, in the case in which p is 

cos / log - j , Poisson's equation is invalid. Instead of cos f log - ) 

we may take, without making any material difference, cos ( log - ) 
plus any constant or differentiate function. If, then, you take a 
solid sphere of density 2 4- cos I log - ) , say, and suppose it divided 
into thin concentric shells, then, by Newton's theorem, the attrac- 
tion due to a shell is zero at a point inside it and therefore 
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da 



= at the centre. The mass of the shells between the centre and 
or 

the point (, r/, ) is 




FIG. C. 



r 

\/2 + cos (log i 



dt. 



Since the attraction at (, 17, ) is the same as if the whole of 
this mass were concentrated at the centre, the attraction is 



and is along the radius vector. Now at the centre 




- , +tao - 

3 \/10 i V b r 

As r is indefinitely diminished, - tends to and cos 1 log - j 

d*u 
makes an infinite number of oscillations. Therefore r-j- is non- 

existent at the centre. Moreover the radius vector may be taken 
in any direction. Hence, at the centre, 

d*u &u d*u 
3? ' dy* ' 3? 

are all non-existent and consequently the equation of Poisson is 
meaningless. 
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For the sake of its historical importance I proceed to give 
you Gauss's proof of Poisson's equation : 

Let /> /(, 17, ) Then the .^-component of the attraction at 
the point (x, y, z) is 

X(x 9 y, 



* f tf-*)/(*,iy,0<ft 

= 1 

J {(-*r + (r/-?y)* + (f- 



where dtf is the element of volume of the sphere and the integral is 
taken over the whole sphere. 
Also 



v 
A (x + e,y,z) 



f ( # e) /(, -/;, ) dt 

j { (x e) 2 " * ( r tj iy)* + (f J 



the integral being taken throughout T 9 the original volume occupied 
by the sphere. Suppose the sphere to be displaced a little to the 
left so that its centre is at P 19 i. e. ( - e, 0, 0). Let T { denote the 
volume occupied by the sphere in its new position. Then, as is 
clear from th^ accompanying figure (7), 




FIG. 7. 



Now, put 
Then 

X(x + e 



=f 



(3) 
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Also 



where Z is the integral in (1) taken over T , A the same integral 
tfaken over 0, l v the integral in (3) taken over 5P and X, the same 
integral taken over r 

The value of an element of volume in 6 and that of an 
element in 0, are respectively e cos dS and - e cos a c?$ where a is 
the angle made by the external normal with the x - axis. 
Therefore in' the limit, 



/v 

taken throughout the surface 8 of T {} ; 
also in the limit 



e * fS <9 



taken throughout T f) . 
Therefore 

0.A O %v t v *- vt> , . i v ~/ Jo 

/)*S ? 

and similarly for 

and 

Qy 

Therefore 

^TT O 

Since 

r^/ =r 2 sin 



the volume integral in the above equality will tend to zero as the 
radius of the sphere is made to tend to zero ; and similarly 

- 1 ^ dtS will approach -4ir/ (#, y, z). 
Thus the equation 

holds. 

I proceed now to consider the generalization of Poisson's 
equation given by Dr. Petrini. Dr. Petrini was the first to note 
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clearly that the continuity of p was not sufficient for the validity of 

cos * 
Poisson's equation. Taking P= - - , he showed * that Poisson's 



equation fails because ^ * s + while r-y and r-j- are bothequftl 
to - OD . Dr. Petrini was thus led to the following generalization : 

lim r 1 fdV (x + h }) y, z) _ ^J^Jj^j \ ff)F (x,y + li^ z) 
= L ^ L 3F ' dx * + h, I dy 



<1V (x, y, z) 



dz -- 81 

provided that while A p h^ h. A are made to tend to zero in any 

manner. 7-^,7-', r should tend neither to zei'o nor to infinity. Dr. 
' /4 2 h.;, h.j 

Petrini proved the validity of his generalization for /*=r ~ and, 



indeed, for every case where /> is continuous. But he did not con- 
sider the cases in which p has a discontinuity of the second kind. 
Such cases have been studied by me in the paper referred to above. 
I have proved, for example, that when 



g -1 , 

] du 4;r f f. 1 -i l\\ 

- 5T = -- 7= 1 cos I log - + tan ~ I f 
r dr ^/IQ L \ G r 37 J 

which shows that the limit in the generalization of Petrini will 
exist and be equal to the limit of 

-4ir ( /, 1 . - 1 U /, 1 , - 1 U 

- ^r J cos I log 4 tan ~ I + cos I log -- -f tan - I 

V/10 I V ^i 3/ V * 3/ 

* cos (log L + tan' 

if the latter limit exists. Now put aA,=A 2 and fth^h^. Then 

1 ~ l 1 
denoting log r- tan - by log H, the expression within the crook- 

/&| O 

ed brackets becomes 



i For another method than Petrini'g, see rny paper in the Phil. Mag., Vol.. 
34, 1917. 



cos 
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(log H) + cos (log H + log - j + cos (log H + log ~ \ , 
i. e., cos (log H) \ 1 + cos (log - ) + cos ( log - \ j 

- sin (log H) | sin Nog -J + sin (log - J | . 

Now the only way in which this expression can have a limit 
is by the vanishing of both 

j 1 +cos (log -J + cos (log j > and j sin (log- W sin (log -H 
which is possible only when a and ft have respectively the forms 



m and n being integers. Thus, it is proved that the generaliza- 
tion of Petrini does not hold for every discontinuous distribution. 

Another generalization of Poisson 's equation was given by 
Professor C. W. Oseen of Upsala. As I have already taken a 
good deal of your time this evening, I will take up the consider- 
ation of that generalization to-morrow. I will conclude the 
lecture by discussing briefly the equation of Coulomb and Poisson 
and Helmholtz's equation. 

I have studied the* equation of Coulomb and Poisson in three 
papers one of which appeared in the Bendiconti del Circolo Mate- 
matico di Palermo in 1917, another in the Philosophical Magazine 
in 1918 and the third in the Bulletin of the Calcutta Mathematical 
Society in 1919. If round P, any point on the surface, a small disc 
be described lying on the surface, then for investigating the validity 
of the equation , ^ ~ v 

_. + _ = - 47H7 , 

an, dn 9 

we need only take into consideration the part of V due to the disc. 
Among other results, I have proved that, denoting by r the 
distance of any point of the disc from P, if 

(T = cos x(r) when \(r) <^0 
and, further, if 

X (r) 1 log i 
r 

dV dV 

then ^ , are both non-existent and the equation of Coulomb 
on i dn z ^ 

and Poisson fails. 

Results similar to those mentioned above have been obtained 
by me for the validity of Helmholtz's equation in a paper which 
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appeared in the Proceedings of the Benares Mathematical Society 
for 1920. I have shown, for example, that, if the strength 



. 1 

(T = COS log - , 



then W + (} and W _ are both non-existent and therefore the equation 
fails. 



FIFTH LECTURE. 

[Delivered on the 4th February, 1921,] 

THE PARTIAL DIFFERENTIAL EQUATIONS OF THE 
THEORY OF GRAVITATIONAL ATTRACTIONS, AND 
THOSE OF ELECTROSTATICS AND MAGNETO- 
STATICS (Continued). 

Yesterday I promised to lay before you the generalization of 
Poisson's equation given by Professor C. W. Oseen. The general- 
ization ] is fairJy simple and will prove as useful as Poisson's 
equation. It is this : 



s n^ 

where dv is an element of the volume o s , ds an element of the 

surface 8 which bounds the volume H K and -^- is the differential 

on 

coefficient with respect to the normal taken outwards. The only 
condition imposed here on p is that it should be finite and 
integrable. 

It will be remembered that in proving the theorem 

we have to impose more restrictions on p than that it should be 
finite and integrable ; for example, in Gauss's proof, it is supposed 
that p is continuous and also that p K , i> 9 , p z exist and are conti- 
nuous. We see at once how much more simple Professor Oseen 's 
generalization is. 

Starting from the equation, known in Integral Calculus as 
Green's theorem, 



one can deduce the modified theorem of Green, viz., 



1 See Oseen' s paper, " Ueber einen Satz von Green und ueber die Definitionen 

von rot und div " / Rendir.n'n.ti d*1 Oi.rr.nln vn.ntevnnt.ifin di Pnlrrtnn t,. ?lft. 
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where // stands for i>(P f ) and V for F(P / ) } P' being the variable 
point (, v/, c), P the fixed point (x, y, 2), and r the distance PP'. 

This result follows from Green's theorem which holds when 
vH7, r*F are existent and integrable. To show this, i take U to 

* )e pyJ> Now pp becomes infinite at the point P. In order 

to make Green's theorem applicable, we enclose the point P within 
a small sphere u' of radius a. Inside the region bounded by this 
sphere and #, U and ite differential coefficients are finite and 

continuous. Hence, since - satisfies Laplace's equation. Green's 
theorem applied to this region gives 



\ 1 r*F'rf?>- \ - 
1 r j r 






_f {!!. r .4 

~" ] I r dn ()u 






where r/s' is an element of the surface 8' of the sphere (2'. Now 
let a tend to zero, then, dv being of the order r*dr, the second 
term in the left-hand side of the above equation is negligible. 
Similarly, as ds f is of the same order as r*, the third term in the 
right-hand side is negligible. Also v*r'=-47iy/ according to 
Poisson's equation. Hence we have in the place of (1), 



. . . 

which becomes, in the limit when a tends to zero, 






This result is the modified theorem of Green and it says that 
the value of F at any point is the potential due to (1) a distribu- 
tion of volume density P', (2) a certain distribution on the bound- 
ing surface and (3) a certain double layer on that surface. Thus 
we see that if there is a solid with a known volume distribution 

and prescribed values of F and ^ on the bounding surface, then 

on 

by applying the modified theorem of Green we can at once put 
down its potential F at any point. 

The question arises : Can we deduce the modified theorem of 
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Green by using Oseen's generalization instead of Poisson's equa- 
tion ? The question lias been answered in the affirmative by Oseen 
who has deduced the modified theorem of Green by applying 
vector analysis and by employing the following theorem : 

If <p and $ are two functions, finite and continuous together 
w}th their first differential coefficients, then 



f l' <ty I' 

\(graci $ x grad ty) dv W r- ds ir\ 



The procedure adopted by Oseen is rather complicated, since 
the second differential coefficients of <l> may be non-existent, but 
the deduction is an epoch-making result. 

The equation of Laplace, which 1 have frequently mentioned 
in some of my previous lectures is intimately associated with 
problems in Electrostatics and Magnetostatics. Thus we come 
across this equation when considering the case of a solenoidal 
distribution of electricity and the most important problem associa- 
ted with the equation is Dirichlet's problem, viz. to find the 
potential at points inside a region when the value of the potential 
is given on the boundary. 

With reference to Dirichlet's problem, two theorems have 
been formulated and proved under certain general conditions. 
These theorems are : 

(1) (The Theorem of Uniqueness) : The solution of Diriehlet's 
problem is unique. 

(2) (The Existence Theorem) : The prescribed value on the 
boundary is arbitrary and corresponding to each value there exists 
a solution. 

That both the theorems hold, was at first guessed from certain 
physical consideration*. For example, in the case of a conductor 
connected with the earth and influenced by external electric 
masses, the value of the potential of the induced distribution on 
the surface is prescribed, being equal in magnitude, and opposite 
in sign, to that due to the masses, and the uniqueness and 
existence of the potential is obvious to a physicist whatever the 
masses may be. But attempts were made to base the two 
theorems on something more rigorous. 

Diriohlet formulated the principle, known as Dirichlet's 
principle, that, of all the functions, F, which inside the region 
bounded by a closed surface remain continuous together with 
their first differential coefficients and take a prescribed value on 
the surface, there must exist one (or more) for which the integral 
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taken over the whole region has a minimum. Accepting this 
principle, we get from the calculus of variations the result that V 
must satisfy Laplace's equation. 

But Weierstrass pointed out in 1860 that Dirichlet's principle 
was not always valid and that all that could be said was that tjie 
integral has a lower limit and that it may not have a minimuvn. 
In 1899, Hilbert established the validity of Dirichlet's principle 
on the supposition that the prescribed boundary value is con- 
tinuous together with the differential coefficients of the first two 
orders. 

Assuming that in a given case the solution of Dirichlet's 
problem exists and is unique, the question arises of actually 
obtaining this solution. I propose to devote the remaining time 
to-day to the discussion of a number of methods. These are (1) 
the use of Green's function, (2) the use of fundamental functions 
based on a specially suited co-ordinate system, (3) the method? du 
balayage, (4) the method of arithmetic means of Poincar6 and (5) 
the use of integral equations. , 

Adopting the procedure used in deducing the modified theorem 

of Green from Green's theorem but taking U to be G + - instead 
of - , we obtain for V(P) the simple expression 



-ij'i 



here G is regular everywhere inside the region, where it satisfies 
Laplace's equation, and equals - - at every point on the surface ; 

G H is called Green's function for the particular region. The 

problem of obtaining the solution of Dirichlet's problem has been 
thus reduced to that of obtaining Green's function which may 
be interpreted as the potential due to a unit of electricity at P 
plus the potential due to the distribution induced on 8 by this 
unit when 8 is connected with the earth. In some cases Green's 
function is easily obtained by applying Lord Kelvin's theory of 
images. For example, let 8 be a spherical surface of radius 
a. Then the value of Green's function at any point Q is 

a 1 J. 
~ R' P'Q*PQ* 

where R denotes the distance of P from the centre of the sphere, 
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FlGL 8. 

P' the inverse of P with respect to the sphere; in fact = units at 

K 

P' form the image of a unit at P. 

I proceed now to consider at some length the second method, 
viz. that in which use is made of fundamental functions based on 
a specially selected system of co-ordinates suitable for the parti- 
cular boundary surface. 

As an example of this method, consider the case in which the 
bounding surface is a sphere. Then, transforming v*F = into 
the polar co-ordinates r, 0, <p, r = a being the bounding spherical 
surface, we have 



fl / 

r>V 



1 fl / . fiV I 



We start by obtaining a typical solution of this equation in 
the form of a product of three factors, one a function of r alone, 
another a function of alone and the third a function of <J> alone. 
Such a solution which is the desired fundamental function, is easily 
found to be 

cos mp 
r"/V" (cos 0) or 

sin m<p 

where 



The solution of Dirichlet's problem is given by 
F= 2 2 rJ> n (A Htm cos m$ + B,, t m sin 

w = 

where the A 's and B's are constants which have to be determined 
by the condition that on the boundary, r=a, V takes the prescrib- 
ed value, say /(0, ). Now in actual practice the determination of 
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those constants is not an easy affair. Gauss had given ' a method 
for expanding any arbitrary polynomial (in the co-ordinates of a 
point on the spherical surface) in a series of spherical harmonics. 
Tf ty(x, y, z) be the polynomial of degree n in .r, y, z where 

a 2 + ?/* + z 1 = n*, 
then Gauss put 

$(x> y, c) = N B -f r* #,-., + r* S n . t + ..... 

S m standing for a solid spherical harmonic of degree m, and, by 
operating on the above equation with r 2 , r* 3 etc., obtained a 
number of simultaneous equations by solving which the $'s were 
obtained. But no compact and general formula for writing down 
the expansion had been obtained until about twenty years ago I 
gave it in a paper published in the M esw.nger of Mathematics, Vol. 
30, in 1900. 



The result is 2n _ L > ul] 



f 



, n 

U 

J r \ 



rv ' *' " ; ~ v ' M 2,1 (2^-2) . . . 2 . (2/4-2,* 4- 1) ('in- 2/i- 1) ...3.1 
where 

tt =- fl = A = /? v 5 ^ + + 

The above-mentioned result was obtained by me at first by a 
method different from that of Gauss whose method was applied by 
me to get the result in a paper published in the Mathematische 
Annalen, Vol. 75, in 1912. 

If / (0, <t>) is expressed as a function of x, y, z, the co-ordinates 
of a point on the boundary, that function can, in its turn? be 
expressed as a sum of polynomials each of which may be expanded 
by my formula, and thus an expansion for /(#, 0) can be obtained. 
With the help of this expansion, the solution of Dirichlet's problem 
can be put down at once as was done by me in my paper of 1900. 

As a second example of the method of fundamental functions, 
consider the case in which the boundary is a cylinder. Then, 
transforming 7*F = into cylindrical co-ordinates />, z, $, p~a 
being the boundary, we get 



The typical solution of the above is 

Jcz cos m f> 

e J m (kp) or ; 
sin m0 

I See his *' Gesamrnelte Werke", Bd. 5, p. 380. 
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this is the desired fundamental function, and the solution of 
Dirichlet's problem is given by 

(v / k n z -k n z \ / \ 

A m , M cos m(f> + B n> )M sin mb) I C tlt m e + /), m e ]J m I k n pl 

where the constants, the A'&, $X 6"s and /)'s are determined by 
the boundary condition. 

As the third example of the method of fundamental functions 
consider the case in which the boundary is an ellipsoid. Then, 
transforming r 2 F = Ointo the ellipsoidal co-ordinates A, /<, r, A = 
being the given ellipsoidal boundary 



X L y 1 



we get 



\vhere 



We suppose, as usual, that a typical solution of the above 
equation can be obtained as a product of three factors, each being 
a function of one co-ordinate alone. Thus we get as the funda- 
mental function 

E(\) E(c) E(i'). 
where 



(I \* 
-~\ E(0)= 



which is Lame's equation and can be solved in scries The 
required solution of Dirichlet's problem is given by 

V = 2A m . n E mt>t (A) E MtH M /?,, (r), 

the /Ts being constants which are determined by the condition 
that, for A = 0, V has the prescribed value, say F(/s ') Many 
attempts have been made to obtain a compact and general formula 
for expanding F(n, ') in ellipsoidal harmonics similar to that given 
by me for expansion in spherical harmonics, but so far the attempts 
have proved to be unsuccessful. 

The question of the determination of the most general system 
of three orthogonal surfaces, namely, that of the confocal cyclides, 
for which the fundamental function is the product of three factors, 
as in the cases considered above, was investigated by Professor 
Klein and his pupil, the late Professor Maxime Boeher. (See 
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Bocher's book, " ReihoneDtwicklungen der Potentialtheorie ", 
Leipzig, 1894.) For the case of an arbitrary bounding surface, 
Poincar6 established in 1896 the existence of fundamental func- 
tions, which, however, are not necessarily expressible as the product 
of three factors, each factor being a function of only one co- 
ordinate. (See his paper in the A eta Mathematica. Vol. 20.) ' 

Unlike the previous method, the methode du balayage of 
Poincare ! uses approximation functions each of which has the 
prescribed value on the bounding surface and each satisfies 
Laplace's equation more approximately than the one preceding it. 
Briefly, the method may be outlined as follows : 

Let K be a sphere which encloses the surface # and the region 
UK bounded by it and let the prescribed value /(#, y, z) be such 
that / together with its differential coefficients of the first two 
orders is continuous in K : also assume that r*fU in K, and let tf> 
denote . 



Then obviously v*(0>-/) = in U 5 . Now construct an enu- 
merable aggregate of spheres K (} , K { , K. r . . . such that each sphere 
lies in Q s and fills it in such a manner that each point of O s is 
inside at least one sphere and the spheres become indefinitely 
large in number only in the neighbourhood of ft. Then, if the 
sequence of functions U (} , U^ U. 2 , ... is determined by the condi- 
tions that (I) V 2 / = in K n and [ 7 equals </> in K - K l)9 (2) 
^ a ?7, = in K } and ?/, equals E7 in A r - A r ,, (3) r 2 ^, = in K., and 
t/ 2 equals t/, in K-K^ etc., f/ ^ U { ^ (J. 2 . . . and the sequence 
t/ , C/p [7 2 , . . . converges uniformly to a regular potential-function 
in O, v : the required solution of Dirichlet's problem is 



The case in which r*/ >0 can be dealt with by making it 
dependent on the case considered above by suitable superposition. 

The method of arithmetic means as given by Neumann * 
postulates the convexity of 8 and determines the solution of 
Dirichlet's problem as the potential of a double layer. If 



W ~~2* I 1 ' f>n\r 



1 See his paper. *Sur les equations aux cterivees partielles de la Physique 
Math6matique " (American Journal of Mathematics, Vol. 12, 1890). 

2 See his booklet, "Untersuchungen iiber das logarithraische and New- 
tonflche Potential, " Leipzig, 1877. 
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' etc - 



8 

then it has been proved by Neumann that the series 



s 

converges in il s , satisfies Laplace's equation and, with the approach 
of the point P, at which the potential is taken, to the boundary, 
tends to /. Here C is a constant and equals the value of 
Jim W tt on the surface. 

n-oo 

The method in which the theory of integral equations is used 
is allied to the above method inasmuch as the solution of Dirichlet's 
problejn is obtained as the potential due to a double layer. Let 
the strength of the layer be />, then it is known that 

W e - W>,= -4^(7'), and 



, = 2 



8 

where W, or W\ is the limit to which the potential W of the 
double layer tends according as the point P on the surface is 
approached from the outside or from the inside, P, M, $ having the 
same meanings as in the figure on p. 6. Adding the two equations, 
we have 

2W e = -47r/> (P) + 2 i P (M) ^r- ds . 
s 

But W, is equal to the prescribed value f(P). Therefore P is 
given by the integral equation 



SIXTH LECTURE. 

[Delivered on the 5tli February, ID21.J 

THE PARTIAL DIFFERENTIAL EQUATIONS OF HYDRO- 
DYNAMICS AND THOSE OF ELECTRODYNAMICS 
AND THE THEORY OF ELECTRONS. 

1 propose to speak to-day first on the partial differential 
equations of Hydrodynamies. Less than three hundred years 
ago, Hydrodynamics was not a science at ail. In fact the first 
theorem ever discovered in Hydrodynamics is what is known as 
Torrieelli's theorem. This was first stated in the year 1643. 
Torricelli was one of the greatest engineers of his day. He /ound 
from observation that the velocity of water issuing from a small 
hole in a vessel of water is the same as that of a particle falling 
from the level of the water in the vessel, i.e., # 2 = 2</z, where q is the 
velocity of water flowing out and z is the depth of the hole below 
the water level. The first attempt at a mathematical proof of 
this, was made by Newton. He pointed out that two important 
points were to be kept in mind. (I) The fluid particles in any 
horizontal section of the jet will always remain in a horizontal 
section. (2) As the section increases in size, the velocity in the 
section becomes less in the inverse ratio. This was a crude form of 
the equation of continuity. From considerations like these Newton 
satisfied himself that the theorem of Torricelli was correct. He 
published his arguments in the Principia in 1687. But the 
theorem was first proved completely by D. Bernoulli in 1738. 
In 1755 Eulergave the equations of motion 

Da ^ I d/> Dv ^ 1 dp Dw __ I dp 

Dt p <9.r ; Dt p dy ' Dt p dz 

, 1) d d <1 d 

where = + ^ --+?> +w 

Dt dt dx dy dz 

u, v t w being the components of velocity and p the pressure at 
(x, y, z). The equation of continuity was also given by him. 
It was first published in the Memoirs of the Academy of Sciences 
of St. Petersburg in 1770. It is 

( ] ^(pu) d(pv) d(pw) ___ 
dl fix dy dz 
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With these equations Euler proved Torricelli's theorem in the 
following mariner : 

When the motion is steady, we have - --- = ~- = = . Also, 

dt at dt 

because the only external force is gravity, we have X=Y = 0, 
and Z~g. The equations of motion become 

du du du I dp dv dv dv I dv 

^ -M, + w - . u 4. v + w = _ . ' . 

ox dy dz p dx r/,r dy dz p dy 
Sw dw dw I dp 

u F* + v "T~ f w ^r = r/ ~ - ;r - 

dx dy r/2: p dz 

Since p is constant for ordinary fluids, the equation of con- 

,. .. , du, dv dw 

tinmty becomes +..+ =o. 
dx dy dz 

Further, since u, v, w are the component velocities in an irro- 
tatiowal motion, we have 

dd> d</) d(f) 

U = i- ?' = ^ 1/5 L 

^ r\ >' t w \ 
dx dy dz 

,. da d*<I> dv AII i 

Hence -- = ; -- = ~- . ATIQ we have two similar equations. 

dy dydx dx ^ 

du flr Ow; 1 a ^ f d / _ 1 5;? 

.*. u +v + w = --- - Or i (u' 2 } = - . 
dx dx dx pdx ' *dx ' pdx 

Similarly, j 1 W )aB - J ^ , and J 1 (g*) = ? - * ^. 
(9t/ p <9// ^ ^2 p dz 

Multiply these equations by dx, dy, dz respectively and inte- 
grate; then 



which is a relation between the pressure and velocity. 

Since the hole is small, we can neglect q* at the level of the 
water. This gives 



P 

.-. J j = ^llf 4. ^. Rut ll2! o, at the hole ; 

^ P 

.'. at the hole, | g*gr, which proves the theorem. 

It will be noticed that the equation of continuity, viz., 

du dv dw> K r x A - i #*< d*< 3V 

* == become8j for irrotatlo ^ motion, -- + -2 , -. o, 
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which is the same as Laplace's equation about which I spoke 
yesterday. Problems in Hydrodynamics can therefore be treated 
in the same way as problems in Electrostatics or Heat. 

For rotational or vortex motion, the equations of motion 
remain the same as before. We have also 

rt> dw dv du dw nt , dv du 
2 = T--~3-; 2^ = --- ; 2 = ~-- ; 
dy dz dz fix dx ay 

where , ^, are the components of tho curl. 

Then 

j^/\ <^ *)<]u du 
Dt \p/ p dv p dy p dz 



1)1 \ p / ~~ p dx p dy p dz 

We can solve these equations in terms of a scalar potential <P 
(no longer a velocity potential). 

CCC ft 1 

If 0=-\\\ f- dx'dy'dz', F=\\\ dx' dt/' dz' > etc., 

* 



r being </(x-x')* +( 

, . du dv dw 

and (9= -T- + V- --^- ; 

d o?/ dz 

then 

a* s# aff 



_ d0 dF dH 
V " By 4 d^""8F' 

~ ? ?~ 

82; (9x dy 

We can thus express u, v, w in terms of #, jP, (?, .?. The anal- 
ysis is very much the same as that which is used for problems on 
the motion of electricity. This is an important point of view 
and is of special importance in the treatment of viscous fluids. 

The first attempt to study the motion of viscous fluids was 
made by Navier in 1827 and the investigation was extended by 
Poisson in 1831. But all the investigations which treated the 
problem from the point of view of molecular theory, took the mole- 
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cules as mere points, although they were supposed to be of finite 
size at the start. These investigations were therefore only fruitful 
in creating interest, but Iwere not of much scientific importance. 
The first person to dea with the problem with any rigour was 
De St. Venant, who is famous for his work in the theory of 
Elasticity. The equations of motion are 

Du 1 dp v de , 



. 

IH p dy 3 dy 

Dw 1 dp v d& 2 

= Z - - -L + - + v? V 
J)t p dz ,5 (h 

where & is the dilatation and equals ~~ + . + . 

(9x dy dz 

These equations have a meaning only when the first and the 
secoijd differential coefficients of n, v, ir exist and are finite. The 
continuity of these had also been postulated, in deriving the equa- 
tions. Prof. Oseen showed that the motion can be determined even 

in the case when U? ~ , ' , etc.. do not exist at all. He 

dx ' dt dx 9 

critically examined the equations, just as I have done for the" 
conduction of heat l and obtained 



as the equation of the phenomenon of linear conduction. 

Oseen 2 has obtained certain very interesting results. I read 
out to you a translation of his own views on the subject : "The 
differential equations of Navier admit of being transformed into 
integral equations in which the differential coefficients 

dy dv dw &u Wu dp dp dp 

W ' dt ' ~dt ' 5? ' % a ' 6tClJ dlr,'fy 9 dz' 

do not occur. Indeed it can be shown that the integral equa- 
tions possess systems of solutions for which these differential 
coefficients do not at all exist. One must then ask whether these 



1 Constitution of Matter and Analytical Theories of Heat " (Abhandbungen 
der Gesellficha/t der Wissenschaften zu Gottingen. Bd. 2, 1903). 

2 " Ueber die Bedeutung der Jntegralgleirbungen in der Theorie der Bewe- 
gung einer reibenden, unzusammendriickbaren Flussigkeit " (Arkiv for Malematik, 
Astronorni och Fynilc, Bd. 0, 1010). 
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solutions of the integral equations, which are not at the same time 
solutions of the differential equations, possess only a mathematical 
importance or whether some among them have a physical import- 
ance. I undertake the task to show that the latter is the case and 
that therefore not the differential equations of Navier but the 
integral equations should be considered as expressing adequately 
the physical hypotheses which lie at the basis of the theory of the 
motion of a viscous fluid. 

"The direct derivation of the integral equations of Hydro- 
dynamics without the intermediation of Navier 's differential equa- 
tions. appears to me to be of interest from two points of view. 
First, obviously the theory gains thereby in generality, secondly, 
the theory gains much in simplicity. In order to see this, we need 
only consider how complicated the considerations are which will be 
necessary to decide under what circumstances the potential func- 
tions will possess differential coefficients of the second order. 
These questions and others of the same kind occur in Hydro- 
dynamics if one starts from Navier's differential equations. If one, 
however, obtains directly the integral equations of Hydrodynamics 
all these questions lose their importance. A similar simplifica- 
tion is obtainable in a similar manner in some branches of mathe- 
matical Physics. I may mention as examples the theory of poten- 
tial, and the theory of the conduction of Heat." 

I give you now the integral equations of Oseen. They bear a 
resemblance to the equations for the vibrating string with an 
infinite number of edges. 

From the principle : 

Amount of motion coming into a closed volume 

= amount of motion due to extraneous forces, 
we get 



,. 

\(/ + A/ - , ) fo- \ dr\ X' 



ete + drX 9 (r) <fe = 



and two similar equations ; 

i TT/ *r I fa & u d u \ 

where X'^X- [u +v -~ + w } , 

V dx dy <1z / 

and X n = Xa, cos nx f X y cos ny + X a cos nz. 

Here X a , X y X, are stresses acting in a plane perpendicular 
to the oj-axis. In the case of a perfect fluid the stress is parallel 
to the axis. In that case X y and X, are zero. Thus in the case 
of viscous fluid there are nine component stresses to consider. 

The differential equation has thus been replaced by an inte- 
gral equation. There is no great difficulty in the solution of such 
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an equation and we can show that the old solution satisfies this 
equation. 

T now come to Electrodynamics. The equations are 



----, 

dy dz 
da dy 

~--~ , 

dz Or 

dp da 
, r - 5 -, 
c)x dy 

da dZ d Y 

dt * dji ~ ~dz ' 



dc_dY dX 
dt ~ dx dy 

where ?/, v, ?/> are the components of the total current strength at 
(x, y, 2), a, ft, y are the components of the magnetic force, a, b, c 
are the components of magnetic induction and X, Y, Z are the 
components of the electric force. The equations can be solved 
by the use of vectors (F, G, H) and the potential (cf. rotational 
motion in Hydrodynamics) . 

The equations of the theory of electrons are more general, 
and include the equations given above as a particular case. 
Prof. H. A. Lorentz has studied the motion of electrons. The 
equations of motion are 

1 <9X u_dy dp 

c dT P c~dy~~Wz ' 

I dY v_da_dy 
c dt c~~ dz dx * 



IdZ w d& da 

-- + p ~ L. ^ __ 

c dt c dx dy ' 

dX dj dZ_ 
te+ dy*te~ pt 

da dZ dY , 

-^7=-^ 3- , etc., as before, 
dt dy dz 

where p is called the volume density of the electricity and u, v, w 
are regarded as the component velocities of the element of electricity 
which is at the point (x t y, z) at time t. 
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The method of solution is as follows : 

dH dG 

Put a = -r -- ;r-, 

dy dz 

fi _M_M 

P~Tz fa ' 

_{)G OF 
y ~Bx,~Ji/' 

Y= _iaF_3F 

c (tl dx ' 

r=- } i- 

c ftt ~dy ' 

li)U BV 

/J ~ --- T?~ V 

< dt dz 

where, F, G, H , are solutions of 

, 1 <VF u 



, . 1 d// w ^ 

and v 2 //- - - + /> -=0 

r* r)r r 



als 



, dF dG dH 1 (9T 7 n 

and ^ ^ ^_ =0. 



r 



There is no difficulty now in getting a solution. If /> = 0, we 
have the equation of wave propagation. Prom this solution as a 
basis we can build up the complete solution. 

Partial differential equations are found to be quite serviceable 
for all ordinary purposes. But in a rigourous treatment of any 
problem they have to be relegated to a secondary place. The 
work of the last 120 years is of importance as it shows that 
most results and processes dealing with differential coefficients are 
valid only under certain conditions. Thus Taylor's theorem is not 

1 

always valid. For example, the expansion of e ^~ a ^ according 
to this theorem, viz.. 

r^-i 2) <*->", 

n ' 
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is wrong because every member of the expansion is zero. This 
discovery, due to Cauchy in 1815, astounded Lagrange who believed 
that if the differential coefficients of all orders existed and were 
finite, the function would be expansible by Taylor's theorem. But 
1 have mentioned to you an example which is in direct contradic- 
tion to his belief. However, in spite of its limitations, the partial 
differential equation will continue to be useful. 

Now I will briefly recapitulate the essence of what I have said 
in criticizing the common use of partial differential equations. 
If we have a string with an edge, the first and second differential 
coefficients do not exist at the edge. The equation of motion be- 
comes meaningless in this case. If we have the initial state of heat 

/(.r) = r for x > 0, and f(x) = - .r for x < 0, 

the first and the second differential coefficients are non-existent at 
the origin. Consequently, the equation of the linear conduction 

of Heat -r- has no meaning for f = 0, ,r 0. In Electrostatics. 
w d.v 2 

Magnetostaties or (Gravitational Potential, if /o = cos log-, r*K 

does not exist. In Electrodynamics too a similar remark applies. 
In Hydrodynamics, as Prof. Oseen says, it is not always the 
differential equations but the integral equations which truly express 
the fundamental physical hypotheses. 

Ifc is quite possible that at some time in the future we shall be 
so advanced in oui thoughts that a differential equation will appear 
to us as but a very crude mathematical instrument, and we shall 
discard it in favoui of the more powerful and more refined integral 
equation. 

It has been a real pleasure 1 to me, gentlemen, to come and 
deliver lectures before you, and 1 shall be much pleased if 1 find 
that my lectures have been of use in fostering a spirit of research 
here. 



